Gradient and Amplitude Scattering in Surface-Corrugated Waveguides 
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We investigate the interplay between amplitude and square-gradient scattering from the rough 
surfaces in multi-mode waveguides (conducting quantum wires). The main result is that for any 
(even small in height) roughness the square-gradient terms in the expression for the wave scattering 
length (electron mean free path) are dominant, provided the correlation length of the surface disorder 
is small enough. This important effect is missed in existing studies of the surface scattering. 



1. The problem of wave transport (both classical and 
quantum) through the guiding surface-disordered sys- 
tems has a quite long history and remains a hot topic 
till now (see, e.g., Refs. 1 2 3 4 5 6 7 8 llfillTlllTlTlll 
EsllT^IT^ and references therein) . One of the main tools to 
treat this problem is the reduction of the wave/electron 
surface scattering to the bulk one in such a way that the 
latter can be described by an effective Hamiltonian with 
a complicated potential, however, with flat boundaries. 
Applying this approach, one can reasonably discriminate 
between the so-called amplitude and gradient scattering, 
and analyze their interplay explicitly. 

To the best of our knowledge, originally the idea of 
this approach was discussed by Migdal^. After, it was 
frequently used in the theories of classical and qua ntum 
wave/electron scattering, see, e.g., Refs. |q|7 p8l9lll0llllllia . 
IT^ But in the majority of themSiSiSiiS the study was 
restricted to the lowest order in the roughness height a. 
Other methodsiiiiSiiiiiiiSiiSiii were mainly based on the 
principal assumption that the surface roughness is suffi- 
ciently smooth. 

In this contribution we present the theory of wave scat- 
tering from rough surfaces, which takes into account both 
the amplitude and gradient scattering mechanisms. The 
important point is that we do not assume any special 
restrictions on the model parameters except for general 
conditions of weak scattering. The latter provides us with 
an appropriate perturbative approach in scattering po- 
tential, however, is not restricted by the smoothness of 
surfaces. 

The main attention is paid to the competition between 
the amplitude and gradient scattering. One of our main 
results is that at fixed r.m.s. roughness height a, the less 
correlation length Rc of a random surface profile, the 
larger contribution of the gradient mechanism. Thus, by 
passing from the smooth to white-noise profiles, the gra- 
dient scattering becomes to prevail. We have analyzed 
this crossover and obtained the estimates describing the 
transition to the dominating gradient scattering. In par- 



ticular, we show that this transition is located within the 
region of small roughness slopes where a/Rc ^ 1. 

2. In what follows we consider an open plane wavegui- 
de (or conducting quasi-one-dimensional quantum wire) 
of the average width d, stretched along the x-axis. For 
simplicity, one (lower) surface of the waveguide is as- 
sumed to be flat, z = 0, while the other (upper) surface 
has a rough profile z = w{x) = d + a£,{x) with (^(x)) = 
and (f^(a;)) = 1. The average (...) is performed over 
different realizations of a statistically homogeneous and 
isotropic Gaussian random function ^{x). We also as- 
sume that its pair correlator {£,{x)^{x')) = yV{x — x') de- 
creases on a scale Rc, with the normalization W{0) — 1. 
The roughness-height power (RHP) spectrum W(k.j;) is 
defined by 



dx exp{~ikxx) W{x). 



(1) 



Since yV{x) is an even function of x, its Fourier transform 
is even, real and nonnegative function of k^. The 
RHP spectrum has maximum a.t k^ = with 1^(0) ^ Rc, 
and decreases on the scale 

In order to analyze the surface scattering for our 
model, we employ the method of the retarded Green's 
function Q{x,x' ; z, z'). Specifically, we start with the 
Dirichlet boundary-value problem 



g2 q2 



^5{x-x')5{z~z'), (2) 

Q{x, x; z = 0, z') — Q{x,x' ; z = w{x), z') = 0. (3) 

Here the wave number k is equal to cj/c for an electro- 
magnetic wave of the frequency oj and TE polarization, 
propagating through a waveguide with perfectly conduct- 
ing walls. As for an electron quantum wire, k is the Fermi 
wave number within the isotropic Fermi-liquid model. In 
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order to express the surface scattering as a bulk one, we 
perform the transformation to new coordinates, 



•^old 1 



ZoMd/[d + a^ix)], (4) 



in which both waveguide surfaces are flat. Correspond- 
ingly, we introduce the canonically conjugate Green's 
function, Gnew — \/w{x)w{x')Qoid and omit the sub- 
script "new" in what follows. As a result, we arrive at 
the equivalent boundary-value problem governed by the 
equation with a "bulk" perturbation potential. 




w{x) I 
w'^(x) 
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Q{x,x';z,z') 



= 5{x-x')5{z-z'), (5) 

Q{x,x' ■,z = 0, z') = Q{x, x'; z — d, z') — 0. (6) 

Here the prime stands for the derivative over x. 

We emphasize that Eq. Q is exact and valid for any 
form of 'w{x). As one can see, the scattering potential 
depends both on the roughness profile (J^{x) and on its 
gradient a^'(x). Moreover, the potential contains the 
term with the square gradient tT^^'^(a:). This term is 
proportional to tr^ and for this reason was neglected in 
all previous studies of transport properties in the surface- 
disordered waveguides. However, as a matter of fact, the 
square gradient introduces the operator V{x) = £,''^{x) — 
iC'^ix)), which plays a special role. Its pair correlator, 

{V{x)V{x')) = 2{^'{x)e{x')f = 2W"\x - x'), (7) 
determines the square- gradient power (SGP) spectrum 



dxexp{~ik^x)W" {x). (8) 



One should stress that although by integration by parts 
the power spectrum of the roughness gradients <t^'{x) can 
be reduced to the RHP spectrum W{kx), this is not pos- 
sible for the SGP spectrum T{kx). This very fact reflects 
a highly non-trivial role of the square- gradient scattering, 
giving rise to the competition with the well known am- 
plitude scattering, in spite of the seeming smallness of 
the term cr^^'^(a;). 

To proceed, we pass from Eq. ij^l to the Dyson-type 
equation, performing the ensemble averaging with the 
use of the technique developed in Ref. |3 The method al- 
lows one to develop the consistent perturbative approach 
with respect to the scattering potential, which takes ade- 
quately into account the multiple scattering from the cor- 
rugated boundary. After quite cumbersome calculations 



we have obtained the average Green's function which in 
the normal-mode representation has the form 



{Q{x, x'] z, z')) = 

^ exp(ifc„|a; — x'\ 
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\x — X 
2Ln 



(9) 



Here fc„ — ^k"^ — [im/d)'^ corresponds to the unper- 
turbed lengthwise wave number k^, and Nd — [fcd/Tr] 
is the number of propagating modes (or conducting elec- 
tron channels) determined by the integer part [. . .] of the 
ratio kd/n. 

3. Our interest is in the attenuation length or total 
mean free path Ln of the n-th mode. Its inverse value 
is given by the imaginary part of the proper self-energy 
and, in accordance with the form of the scattering poten- 
tial, consists of two terms describing different scattering 
mechanisms, 



1 

-(1) 



(10) 



The first length Ln is determined by the expression 



(nn/d)^ ^ (nn'/d)^ 



(1) 



n' — l 

X [W{kn + fc„') + W{kn - fc„/)] , 



(11) 



Its diagonal term is formed by the amplitude scattering 
while the off-diagonal terms result from the gradient one. 
Eq. Hll|l coincides with that previously obtained by dif- 
ferent methods (see, e.g., Ref. Q). 

The second length Ln is associated solely with the 
square- gradient mechanism due to the operator V(x), 
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(12) 



n' — l nn' 

Its diagonal term controls the intramode scattering, 
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{■Ku/df 
X [T(2fc„) + T(0)] , 



(2) 

The off-diagonal partial length -^„_^„/ describes the in 
termode scattering (from n to n' ^ n channels), 

1 (7rn/d)2 (7rn7d)2 (n2 +n'2)2 
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(13) 
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kn 

X [T(fc„ + kn 



kn' 
+ T{kn 



(14) 



To the best of our knowledge, in the surface-scattering 
problem for multi-mode waveguides the operator V{x) 
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was never taken into account, and, as a result, the square- 

(2) 

gradient attenuation length L„ was missed in previous 
studies. 

Let us analyze the conditions under which Eqs. I|1U|) 
- (tT^ are derived. We stress that the Dyson-type equa- 
tion for the average Green's function was obtained within 
the second-order approximation in the perturbation po- 
tential. This means that the self-energy in this equation 
contains the binary correlator of the surface-scattering 
potential and the unperturbed Green's function. In 
terms of the diagrammatic technique this is similar to 
the "simple vortex" or, the same, Bourret approxima- 
tioni^. Following the ideas discussed in Ref. one can 
show that this approximation is justified when the broad- 
ening 1/2L„ of the quantum wave number kn is much 
less than both the correlation scale and the spacing 
\kn — kn±i\ ~ \dkn/dn\ between neighboring quantum 
wave numbers. The same conditions also arise due to an- 
other approximation which is the use of the unperturbed 
value kn in the expression for the self-energy, instead 
of the perturbed one. Now we take into account that 
\dkn/dn\ ~ A~^, where A„ is the distance between two 
successive reflections of a wave from the rough bound- 
ary inside the n-th channel. As a result, we come to the 
following conditions of a weak surface scattering 

A„ = 2knd/{Tm/d) <C 2i„, Rc < 2i„. (15) 

These inequalities imply that the electron/wave weakly 
attenuates on both the correlation length Rc and the cy- 
cle length A„. 

As one can see, the expressions (|TT|l and ((T^ - 
represent, respectively, basic contributions from princi- 
pally different surface-scattering mechanisms related to 
the amplitude, gradient and square- gradient terms. It 
should be emphasized that the corrections proportional 
to cr'*, originated from higher order approximations in the 
amplitude and gradient terms of the perturbation poten- 
tial, are smaller than the main contribution (|ll|l under 
the conditions (|15|l . Contrary, the square- gradient terms 
give rise to the cr^-terms in Eas. (|13|l - (|14|) which should 
not be neglected due to a specific dependence on the cor- 
relation length Rc- Note that Eq. (|15|l implicitly includes 
the requirement for the surface corrugations be small in 
height {a d), but does not restrict the value a/Rc of 
the roughness slope. 

4. Since Ln'^ and L^n"^ depend on as many as four 
dimensionless parameters (ka)'^, kRc, kd/ir, and n, the 
complete analysis appears to be quite complicated. For 
this reason, below we restrict ourselves by the analysis of 
the interplay between Ln'^ and Ln^ as a function of the 
dimensionless correlation length kRc for Nd ~ kd/n ^ 1. 

As follows from Eq. (|ll|l . the inverse value of the 
amplitude-scattering length typically increases with an 
increase of kRc- Specifically, in the case of the small- 
scale roughness {kRc 1 ^ fcA„) we have l/Ln ^ oc kRc- 
Then, within the intermediate region where 1 <^ kRc ^ 
kAn, the increase of slows down, or can even be 



replaced by the decrease for some values of the model 
parameters. Finally, for large-scale roughness and strong 
correlations (1 < kAn kRc) the value of l/Ln'^ again 
starts to increase linearly with k Re- 
in contrast with 1/Ln\ the inverse square-gradient 

(2) 

scattering length 1 /L„ ' reveals a monotonous decrease as 
the parameter kRc increases. At small {kRc <C 1 ^ fcA„) 
and extremely large (1 ^ fcA„ ^ kRc) values of kRc, 

this decrease obeys the law 1/Ln'' oc {kRc)~^, due to 
T(0) ^ 

/^From this analysis it becomes clear that the curves 
displaying 1/Ln^ and must intersect, and one can 

observe the crossover from the square-gradient to am- 
plitude surface scattering. To the left from the crossing 
point {kRc)0 the square-gradient scattering length pre- 
vails, Ln'^ <C Ln"^ . To its right the main contribution is 
due to the well known amplitude scattering, Ln'' <C in ^■ 
If the crossing point falls onto the interval of the small- 
scale roughness {kRc <C 1), its dependence on the model 
parameters is described by 

{kRc)l ^ {ka)n/^/}~d. (16) 

This estimate shows that the crossing point is smaller for 
smaller values of the dimensionless roughness height ka, 
as well as for smaller mode indices n, or for larger values 
of the parameter kd/n. 

In Fig. n we display the dependence of A„/2L„ as a 
function of kRc assuming the Gaussian binary correlator 
yV(a;) = exp(— a;^/2i?^) for random surface profile £,{x). 
The curves are plotted starting from such values of kRc 

(2) 

for which A„/2L„ — 1, according to the first condition 
of Eq. H15(l . Taking into account the second condition re- 
stricting the maximal value of kRc, we plot every curve 
within the range where Rc < 2Ln\ As one can see, 
all curves have the crossover from the square-gradient to 
amplitude surface scattering. The first (lowest) one with 
{kcr)^ = 10^^ has the crossing point {kRc)Q « 0.2 lo- 
cated within the interval of small-scale roughness, and 
the crossover reveals a small dip centered at {kRc)Q. 
The curve obeys the asymptotic behavior {kRc)^^ to 
the left from {kRc)Q due to the main contribution from 

An/2Ln^. Then the quantity An/2Ln^ becomes dom- 
inating in the sum IjlUI) . therefore, the curve starts to 
rise. Firstly, the linear dependence on kRc on the right 
deep-side (where kRc < 1) is replaced with a smoother 
one (for kRc > 1). Finally, for Rc > An (strong correla- 
tions) the linear dependence restores. 

The crossing points of the second, third and fourth 
curves have the values of the order one. Here the total 
attenuation length L„ within the whole small-scale region 

(2) 

is formed by the square-gradient scattering length L„ . 
In full agreement with Eq. 1161) the presented curves dis- 
play that the smaller parameter (fccr)^ the smaller value 
of the crossing point (fciJc)©- 

Note that for all curves in Fig.l the roughness height 
is small, a/d <^ 1. Furthermore, for the amplitude- 
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FIG. 1: A„/2Ln versus kRc at kd/n — 62.5, n = 31 and 
different (ka)^. 

dominated scattering (to the right from the point (fciJc)© 
where An/2Ln^ mainly contributes), the average corru- 
gation slope is also small for all data, a/Rc 1. The 
roughness slope remains to be small at the crossing points 
too, but increases to their left with the decrease of kRc- 



As a result, to the left from the crossing point where the 
square-gradient term A„/2L„ prevails, the slope reaches 
the values of the order one, or even larger for the first tree 
curves. 

In conclusion, we have discovered the principal impor- 
tance of the square-gradient surface-scattering mecha- 
nism which was never taken into account in the liter- 
ature. We have shown that at any fixed value of the 
roughness height a one can indicate the region of small 
values of the correlation length Rc where the new square- 

(2) 

gradient scattering length L„ predominates over the 
known amplitude scattering length Ln'' {Ln'' ^ ^n^)- 
The predominance occurs in spite of the fact that 1 / Ln^ 

is proportional to while l/L^? is proportional to 
a^. This happens since the two lengths are determined 
by the substantially different roughness-height W(kx) 
and roughness-square-gradient T[kx) power spectra, that 
have vastly different dependencies on R^. It is remark- 
able that the square-gradient mechanism prevails in the 
commonly used region kRc <C 1 of a small-scale bound- 
ary perturbation, where the surface roughness is typically 
described via the white-noise potential. 
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